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ABSTRACT 


The performance of convolutional codes in fading channels typical of the 
planetary entry channel is examined in detail. Here the signal fading is due 
primarily to turbulent atmospheric scattering of the RF signal transmitted 
from an entry probe through a planetary atmosphere. The primary motivation is 
in support of the Pioneer-Venus mission although the results should have wider 
applicability. Short constraint length convolutional codes are considered in 
conjunction with binary phase-shift keyed (BPSK) modulation and Viterbi maximum 
likelihood decoding while for longer constraint length codes we consider 
sequential decoding utilizing both the Fano and Zigangirov-Jelinek (ZJ) algorithms 
Careful consideration is given to the modeling of the channel in terms of a few 
meaningful parameters which can be correlated closely with theoretical propagation 
studies. For short constraint length codes we are primarily interested in the 
bit error probability performance as a function of E^/N q parameterized by the 
fading channel parameters. For longer constraint length codes interest will 
center on the effect of the fading channel parameters on the computational require 
ments of both the Fano and ZJ algorithms. In either case the effects of simple 
block interleaving in combatting the memory of the channel is thoroughly explored. 
The approach is analytic where possible otherwise resort is made to digital 
computer simulation. 
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1.0 Introduction 

The use of convolutional codes in conjunction with coherent binary phase- 
shift keyed (??GK) modulation has proven to be an effective and efficient means 
of obtaining error control on the classical deep space channel. Recent work by 
Heller and Jacobs [1] has discussed the performance of short constraint length 
convolutional codes in conjunction with coherent BPSK modulation and Viterbi 
maximum likelihood decoding (cf. [ 2 ] , [ 3 ] ) while Jacobs [V] has treated the per- 
formance of longer constraint length convolutional codes utilizing sequential 
decoding (cf. [5] , Q$] ) . Again these works have been confined to the classical 
deep space channel where the net effect of the channel is simply the introduction 
of an additive white Gaussian noise (AWGN) component. In future planetary entry 
missions, however, the classical deep space channel provides an inappropriate 
model of the actual propagation environment experienced by an entry probe in 
transmitting an RF signal through a planetary atmosphere. Here it is expected 
that the signal will undergo fading due to a number of causes most notable of 
which appear to be: entry probe dynamics; effects of oscillator instabilities 

which become pronounced at the typically low bit rates associated with planetary 
entry missions; multipath reflections off the planet surface; turbulent atmospheric 
scattering of the RF signal within the planetary atmosphere. Our interest is 
primarily with this latter problem although the approach is purposely quite 
general and can be easily modified to apply to more general fading channel 
conditions. Unfortunately little is known concerning the performance of con- 
volutional codes in fading channel environments typical of the planetary entry 
channel. The present study has been addressed to this problem. The primary 
motivation is in support of the Pioneer-Venus mission although the results should 
have wider applicability. In particular, we consider coherent BPSK modulation 
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in conjunction with Viterbi maximum likelihood decoding for short constraint 
length codes while sequential decoding strategies are employed for the longer 
constraint length codes. In the latter case both the Fano and Zigangirov-Jelinek 
(ZJ) algorithms are considered. Careful consideration is given to the modeling 
of the channel in terms of a few meaningful parameters which can be correlated 
closely with theoretical propagation studies. For short constraint length codes 
we are primarily interested in the bit error probability performance as a function 
of E^/N parameterized by the fading channel parameters. For longer constraint 
length codes interest will center on the effect of the fading channel parameters 
on the computational requirements of both the Fano and ZJ algorithms. In either 
case the effects of simple block interleaving in combatting the memory of the 
channel is thoroughly explored. The approach is analytic where possible otherwise 
resort is made to digital computer simulation. 

2.0 Preliminaries 

Before proceeding to a discussion of binary convolutional code performance 
on fading channels typical of planetary entry missions, let us first consider the 
performance afforded by the use of coherent BPSK modulation on the additive white 
Gaussian noise (AWGN) channel. This will not only serve to introduce the results 
to follow but will allow a check on the accuracy of the various channel/decoder 
simulation programs employed. Consider then the communications system illustrated 
in Figure 2-1. Here the binary data stream is applied to a binary con- 
volutional encoder resulting in a binary sequence (actually + l) at its 

output which is in turn applied to a modulator/transmitter producing at its output 
•the analog channel waveform. 

s(t) cos (oj Q t + 0) I x i p(t - i T g ) 


i 


( 2 . 1 ) 





original 

OF POO R 


p AGE IS 

quality 


Data Communications System Employing Binary Convolutional Coding 
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where 'i> is the carrier frequency, 9 the carrier phase assumed known at the 
o 

receiver, and p(t) a pulse-like waveform of duration T seconds satisfying 

s 

T 

f S p 2 (t)dt = l (2.2) 

'0 

so that E can be interpreted as the signal energy per transmitted channel symbol, 
s 

The quantity E is related to the energy per transmitted information bit E, by 

S D 

E s =RE b < a - 3 > 

where R is the rate of the binary convolutional code. We have been concerned only 

with codes for which R = l/n where n is the number of transmitted channel symbols 

per information bit. In some cases the binary information sequence ^a^ must be 

blocked into segments or frames of length N bits to each of which K-l (K is the 

code constraint length) tail zeros are added which, of course, carry no information. 

Such a procedure is not necessary for Viterbi decoding but is required with 

sequential decoding to facilitate decoder synchronization and/or restart after 

buffer overflow. The effect of the tail zeros is to reduce the energy available 

for information transmission, resulting in a relationship between E g and E^ given by 

E S - *b . 

so that if N» K as in most cases of interest, we have E g fa RE^ as in (2.3). 

• Due to the presence of additive channel noise, the channel output is 

v(t) = s ( t) + n(t) (2.5) 


where n(t) is zero-mean AWGN with double-sided noise spectral density N q /2 watts/Hz. 
We will assume that the receiver utilizes a coherent matched filter matched to the 
pulse p(t). . As a result, the receiver output after suitable normalization is the 
sequence of real numbers 5 r^ 



where 
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i = 1, 


2 ,... 


( 2 . 6 ) 


and, as before, is a sequence of + 1 values determined by the actual infor- 

mation sequence and the specific structure of the binary convolution code employed. 
The sequence ^n/^ , on the other hand, is a sequence of independent and identically 
distributed (i.i.d.) random variables with zero mean and unit variance. The 
decoder processes the sequence W to produce at its output the sequence 

where at. is an estimate of the i^h transmitted information bit a. . Let us 
i i 

consider first the bit error probability performance of this system employing short 
constraint length binary convolutional codes and Viterbi maximum likelihood 
decoding. 


2.1 Viterbi Decoding 

Recently Viterbi £ 3 ] has provided a rather complete analysis of binary con- 
volution codes and their performance in typical communication systems. His work 
has shown the utility of the state diagram and associated generating function 
approach to analyzing such problems . We will assume knowledge of this approach 
in what follows. Of immediate interest is the fact, as Viterbi has shown, that 
a tight upper bound on bit error probability can be obtained in terms of the 
generating function T(D, N) (see Ref. [V] for definitions) associated with a 
particular binary convolutional code. In particular, if d^ is the free distance 
of a 'code of rate R and with associated generating' function T(D, N), then a tight 
upper bound on the resulting bit error probability is provided by the expression 


< erfc ( 




dT(D, N) 

<m 


N=l, D=exp 
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An efficient computer program (cf. ra ) has been written to evaluate this bound 

for arbitrary binary convolutional codes specified only in terms of their code 

connection vectors. Odenvalder [jsj and more recently Larsen 00 have tabulated 

a number of good codes which, among other things, possess the desirable property 

of maximizing the free distance d^. for a given rate and constraint length. In 

particular, Odenvalder lists codes for rates R = l/2 and l/3 and constraint 

lengths K ^ 9, while Larsen has extended this tabulation to R = l/2, l/3 and. 1/4 

for constraint lengths K ^ 24. The upper bound given by (2.7) has been evaluated 

for the Larsen codes for rates R = l/2, l/3 and 1/4 and constraint lengths K 4c 10. 

The results are illustrated in Figures 2-2 through 2-4. The simulated bit error 

probability performance for selected R = l/2 and l/3 codes for various degrees of 

receiver output quantization are illustrated in Figures 2-5 through 2-8. The 

random number generator used in these simulations is one developed and tested 

previously at Lincoln Laboratory Jio} . The path memory of the Viterbi decoder 

was chosen long enough to render negligible the effects of memory truncation 

(approximately 10 constraint lengths in each case). In all cases we have used 

uniform quantization of the matched filter outputs. More specifically, assume 

that there are Q, distinct quantization levels £ q^, q^, ..., where Q is some 

** til 

power of 2. Then if I. represents the upper endpoint of the j quantization 

I 

interval j =1, 2,...,Q we shall assume 

V2 = ° 


. Vj ' L 5 

with the additional assignments 


i = 1 2 — - 1 

j — i-i *->••• > 2 X 


(2.8a) 

(2.8b) 


V + '°° ; 


! q = - CO 


The 


quantization levels are taken as the integers +1, +2,..., +Q/2 


(2.8c) 




Bit Error Probability P& 





Bit Error Probability Pb 



Figure 2-5 

Simulated Bit Error Probability 
Performance for K=3> R=l/2 Code 
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If represents the sequence of quantized matched filter outputs then 

if J - 1 ’ 




'fhe receiver quantization is then 

Q 

completely determined by specification of the quantities i., j = 1, 2 , . . . , - 1. 

We have taken in all our studies the uniform quantization grid such that 

j=2,...,Q/2 (2.9) 

Simulation results have suggested that A values as indicated in Table 2-1 are 
close to optimum and will be utilized in all simulations to be reported here 


I. - I. , = A ; 

3 J-l 



Table 2-1 

Approximately Optimum Uniform 
Quantization Spacing 


The results illustrated in Figures 2-5 through 2-8 are in complete agreement 
with those reported previously by Heller and Jacobs [lj and serves as a useful check 
on the accuracy of the simulation program. 


2.2 Fano Decoder 

An important and useful parameter in characterizing the performance of 

sequential decoding schemes in general and the Fano algorithm in particular is the 

critical rate R (sometimes called the computational cutoff rate R ) which for . 
o comp 

coherent BPSK modulation and equally probable signaling is given by (of. [V ] , 

Chap. 6) 

I _ 

1 + / \J P £r|x=l} p ^ r|x=-l} dr (2.10) 


R o 1 log 2 
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where p { 'I'- - is the a posteriori probability of having observed a receiver 
output value r during any one channel signaling interval given that the corre- 
sponding transmitted symbol was x = i fur i = + 1. For the AWGN channel, the 
receiver outputs are described by (2.6) and the sequel with the result 


i r[ x^- = — - — exo 


- | (r - x 


2E 2 

N ; 
o 


( 2 . 11 ) 


Substituting into (2.10) and performing the indicated integration we obtain 

E 


R o 1 loS 2 


1 + exp < - — 


( 2 . 12 ) 


For quantized receiver outputs, the value of R q must be computed according to 


R = 1 
o 


log. 


1 + 



p i ?= q ol X=1 ] 


(2.13) 


where now p is the a posteriori probability that the quantized 

receiver output assumes the j quantization level q., j = 1, 2, ...,Q during any 
signaling interval given that the corresponding transmitted symbol was x = 1 with 
i = + 1. A plot of R q vs. E s /N q is provided in Figure 2-9 for different amounts 
of receiver quantization. In Table 2-2 we illustrate how closely the value of R q 
for uniform quantization approaches that with optimum quantization for selected 


values of E s /N q . This table is taken from Richer [ll] . 

One of the principal problems encountered in the practical implementation of 
sequential decoding schemes in general and the Fano algorithm in particular is the 
variability of the amount of computation required to decode a message block. In 
particular, if is the total number of computations required to decode a message 


block consisting of N bits, then it can be shown [121-113] that the number of 


computations per decoded information bit C=Cjj/N has a probability distribution 
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which is approximately Pareto (cf. Cramer [l4] , page 2 { i8) of the form 

pjc^xj^AX* 0 (2.14) 

where the approximation holds only for large X. The quantities A and a in (2.l4) 

depend primarily on the channel and, in particular, on the quantity R /r through 

its deoendence on E /N . Here R/R is the code information rate normalized by the 

s' o ' o 

critical rate R and is easily commit ed with the aid of Figure 2-9 as a function 
o 

of E /N . For a > 1, the expected value of the number C of computations per 
decoded information bit exists and is finite while for Ot-£s 1 it does not exist. 
Furthermore, it can be shown that a > 1 for R/R < 1 while ail for R/R ^ 1. It 
is for this reason that R q has been called the critical or computational cutoff 
rate of sequential decoding. Our interest in sequential decoding will be confined 
to a particular code and block length. In particular, we consider a K = 32, 

R = l/2 code described by Massey and Costello [ 15 ] . Although not a systematic 
code, this code has the unique ability of recovering the transmitted data sequence 
by simple modulo-two addition of the hard decisioned received data sequence. In 
all cases we consider a block length of N = 250 bits and Q = 8 level quantization. 

The decoder utilizes threshold spacing Tq = 4 bits. In Figure 2-10 we illustrate 

the empirically derived probability distribution of for several values of 

/ 4 

R/R q determined by simulation on the basis of 10 successive transmissions of the 

* 4 

250 bit message. A computational cutoff of 5 x 10 computations per block has 
been -imposed in these simulations. Table 2-3 lists other pertinent results of 

these simulations. 


The number of computations per decoded information bit C can then be determined 
as C=Cff/Nwith N = 250. 
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Empirical Probability Distribution of Computational Count 
for Fano Decoder in AWGN Channel 
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R / R 0 

Block Error Probability 

* 

Quit Probability 

A 

a 


P E 









0.80 

0.00 

0.00 

4.5x10 s 

3.50 

0.90 

0.00 

0.00 

2.5x10? 

2.82 

0.99 

0.00 

0.002 

1.1x10 s 

1.20 


Based upon an imposed computational cutoff of 5 x 10 computations 
per 250 bit block. 


Table 2-3 

Other Pertinent Results of Fano Decoder 
Simulations in AWGN Channel 

2.3 Zigangirov-Jelinek (Zj) Decoder 

The Fano algorithm treated in the preceding section has been used extensively 
in sequential decoding applications. It has been thoroughly investigated and its 
behavior well documented [l6] - [ 20 ] and generally understood. Recently an alter- 
native sequential decoding algorithm has been proposed independently by 
Zigangirov [21] and Jelinek [22] sometimes referred to as the Zigangirov-Jelinek 
or simply the ZJ algorithm. This algorithm differs from the Fano algorithm in 
that it saves all information on path segments which have been previously examined 
for possible future use. It can be shown J23] that the performance of the two 
algorithms are identical in the sense that,- for all practical purposes, the sets 
of nodes examined by both algorithms are the same as are the paths ultimately 
selected. The significant difference is in the search strategies. With the ZJ 
algorithm, a node may be extended at most once, while the Fano algorithm allows 
repeated visits to a node in the process of decoding a received data sequence. 

There is an obvious tradeoff here between the smaller number of computations of 
the ZJ algorithm and the reduced memory requirements of the Fano algorithm. In 
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light of the potential computational advantage of the ZJ algorithm, and the 
relatively low and continually decreasing cost of memory, the ZJ algorithm is 
the more likely candidate for providing near real-time decoding of long constraint 
length convolutional codes. 

The details of a digital computer program for simulating the performance of 
the ZJ algorithm under a variety of channel conditions has been described 
previously in [24] and need not be repeated at this time. The basic ZJ algorithm 
is' characterized by the use of a memory table or "stack" to store node or 
cumulative path metrics associated with paths previously examined and extended. 
Each table entry consists of the identification of the most recent node on each 
path previously explored together with the cumulative path metric up to that node 
in the code tree. It is to be noted that there is exactly one entry for each 
previously explored path. The table entries are ordered in terms of decreasing 
path metrics so that, in particular, a pointer directed to the top of the stack 
will always be pointing to the node which currently possesses the largest 
cumulative path metric associated with it. The ZJ decoder then consistently 
attempts to extend the path through the node currently identified with the top 
entry in the stack. More specifically, the ZJ algorithm proceeds according to 
the following rules : 

1. ) Initialize by clearing the memory table and creating 

one entry corresponding to the root of the code tree. 

The cumulative path metric associated with this entry 

is set equal to zero. 

2. ) Retrieve the entry with the largest cumulative path 

metric currently residing at the top of the stack. 

If the associated node is at the end of the code tree, 
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decoding is completed - announce the corresponding 
N-bit message as having been transmitted. Otherwise 
continue to step 3- 

3.) Compute the cumulative path metrics associated with 
the two successors of the node identified in step 2. 

b.) Entries are created for these two successor nodes and 
inserted in the stack at appropriate positions while 
the entry for the predecessor node is deleted. 

5.) Go to step 2. 

There are several comments appropriate to this basic ZJ algorithm. First, 

observe that since there are a finite number of nodes in the code tree the 

algorithm will eventually terminate although possibly with an erronous decision. 

Indeed, even when a correct decision is eventually made a large number of 

computations may be required. Second, it is to be noted that a decoding table 
N 

with at least 2 entries is required to avoid overflow although under typical 
noise conditions only a small percentage of this number will actually be used. 
For large IT it becomes unfeasible to provide the required memory so that the 
ZJ algorithm will be forced to operate with a reduced table size. Under such 
conditions it is, of course, possible that momentarily atypical noise conditions 
cam result in the only entry associated with ‘the correct path being purged from 
the stack. As a result, the decoder will have great difficulty in extending 
any of the remaining incorrect paths to the end of the code tree. It is clear, 
however, that given enough time the decoder will eventually make a decision 
although in error. If the table size is to be reduced then, it would appear 
appropriate to terminate decoding with declaration of an erasure if the number 
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* 

C of computations per decoded information bit exceeds a fixed preset limit C m&x , 
or equivalently if the total number of computations exceeds NC mav . This will 
surely reduce the block error probability P_ but at the expense of an increase 

lli 

in the quit probability P^. This is tolerable provided that C max is not so low 

as to result in the decoder quitting an appreciable percentage of the time when 

the correct path is in fact still represented in the decoding table. As part of 

N 

the description of the ZJ decoder then we must specify the table size T^ 2 

and the imposed computational cutoff C mfly . Other than some experimental work by 

Richer Jllj there is little information to aid in this choice. In Figure 2-11 

we illustrate the empirically derived probability distribution of for the ZJ 

decoder for several values of r/R q again using the K = 32, R = l/ 2 Massey 

4 

quick-look code. These results were obtained on the basis of 10 successive 

transmissions of the 250 bit message with an imposed computational cutoff of 
4 

5 x 10 computions as in the case of the Fano decoder. Again the computational 
distribution exhibits Pareto behavior although the computation count for a 
specified r/R q is, as expected, reduced considerably over that for the Fano 
decoder. The table size T in these simulations was chosen to consist of a 
maximum of 5 x 10 entries and is sufficient to insure that the probability of 
purging the -correct path is negligible. Indeed, Figure 2-12 illustrates the 
empirically determined probability distribution of the required table size for 
several values of R/R , again for the Massey K = 32, R = l/2 code. This data 
was obtained by monitoring the maximum depth into the decode table or stack 
achieved by the table entry corresponding to the correct path. Interestingly 
enough this distribution again exhibits Pareto behavior although we have not made 
any great effort to pursue this result. It is clear, however, from Figure 2-12 


Here a computation will be taken to mean a single iteration of the steps 1-5 
describing the basic ZJ algorithm. 
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that only moderate table sizes are required to provide reliable decoding in the 
presence of AWGN. Table 2-4 lists some other pertinent simulation results. 



Based upon an imposed computational cutoff of 5 x 10 computations per 
250 bit block and a decode table consisting of 5 x 103 entries. 


Table 2-4 

Other Pertinent Results of ZJ Decoder 
Simulations in AWGN Channel 

It is to be noted that the computational advantages of the ZJ decoder vis-a-vis 
the Fano decoder is offset somewhat by the more complicated nature of a typical 
ZJ computation. In particular, even for decode tables of modest size, the sort and 

merge (or insert) operation represented by step 4 of the basic ZJ algorithm can be 

1 

extremely time consuming if not performed intelligently. As a result, the com- 
putational advantage of the ZJ algorithm can be severely compromised if careful 
attention is not given to the data file organization and/or search strategy of the 
decode table. We have not investigated this topic to any great extent. It is 
felt,, however, that with appropriate care given to such details the ZJ decoder 
does indeed appear a more attractive alternative to the Fano decoder for providing 
near-real time decoding of long constraint length convolutional codes. 
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3.0 Fading Channel Characterization 

The transmitted channel waveform s(t) for BPSK modulation has been given 
previously by (2.2) and the sequel. It will be more convenient in what follows 
to consistently use complex signal representations. In particular, s(t) can be 
expressed in the form 


s(t) = Re 




u (t-iT ) 
o s' 



(3.1) 


where u q ( t) is the complex envelope of the transmitted channel signal during any 
one signaling interval. More specifically, in this case 

u Q (t) = p(t) e^ 0 (3.2) 

where p(t) and 9 have been defined previously in conjunction with (2.1). From 
(2.2) it follows that 

J TS | u o (t) | 2 dt = 1 (3.3) 

0 

The transmitted signal s(t) can be written' in the still more convenient form 

c 1 7 

s(t) =Re ] u(t) e j (3.4) 


where now 

u(t) = \/2 eJ x i u 0 ( t " iT s ) (3*5) 

i 

The received signal is similarly of the form 

C ■ t ■) 

u(t) = Re ^w(t) e j (3-6) 

where now w(t)- is the complex envelope of the received waveform and will be assumed 
of the form 

w(t) = z(t) + n(t) 


(3.7) 
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Here the quantity z(t) is the complex amplitude of the received scattered signal 
component and n(t) is a signal independent zero-mean stationary complex Gaussian 
random process with 

E ^n(t) n(t Q )^ = 0 (3.8a) 

and 

E ^n(t) n*(t Q )^ = N q £(t-t Q ) (3-8b) 

where £(•) is the Dirac delta function and N q /2 is the double-sided noise spectral 
density of the inphase and quadrature components in watts/Hz. We have assumed that 
z(t) can be expressed in the form (cf. { 25 } ) 

z(t) » [~T + a(t)J u(t) (3-9) 

Here -T = X is a complex quantity whose amplitude X is a fixed deterministic 
quantity while the phase tp is uniformly distributed over C-it, jr"] . The quantity 
a(t) is a complex zero-mean Gaussian random process completely described in terms 
of a channel scattering function as described previously in [25] . We will assume 


that is the common variance of the inphase and quadrature components of a(t) 

2 2 r 

(actually = 0^ /2 in the notation of [25j ). Under these conditions, the 

probability density function (p.d.f.) of 

r(t) = | r + a(t) j 

is given by the Rayleigh-Rice distribution 


(3.10) 


V r> ■ 53 exp < - 

= 0 


2cr 2 


0 (j-2. 


r ^ 0 
r <0 


(3.11) 


Hopefully it will not cause confusion that we are using r(t) to represent the 
envelope of the received signal component and to represent the sequence 

of matched filter outputs. 
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where I (x) is the modified Bessel function of the first kind of order zero. If 
o 

^ = 0 (or equivalently f = 0), this reduces to the Rayleigh distribution. 



r >> 0 
r < 0 


(3.12) 


The mean-square value of the Rayleigh-Rice distribution is easily seen to be given 
■by 

E = 2 C' 2 + X 2 (3-13) 

= <r 2 [2 + a 2 ] 

A 2 2 

where a = 'i/ff' allowing ^ to be measured in units of . 

It should be noted, that a number of theoretical propagation studies concerned 
with RF signal fading in the Venusian atmosphere (cf . [2 6] , [27] ) have been con- 
cluded indicating that instead of (3. 9) > the received signal component is described 
by 

z(t) ='a(t) u(t) 

where the amplitude r(t) = | a(t) | is lognormal, i.e., 

r(t) = e %(t) (3.14) 


2 

with %(t) a real zero-mean Gaussian random process with variance and the 

phase 0(t) = arg ^ a(t)^ is likewise a real . zero-mean Gaussian random process with 
variance . It follows that Under the lognormal assumption r(t) has mean 





E ^r^ = exp ^ 


(3.15) 

and mean-square value 


■ 

E ^r 2 ^ = ex P 

( 2 

(3.16) 
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Table 3-1 summarizes some pertinent information on the log-amplitude process *%.(t) 

and the phase process 0(t) taken directly from Q?t| • 

The results of these previously referenced propagation studies have been 

discussed in [ 28 ] where arguments have been advanced in support of the Gaussian 

model represented by (3-9) and the sequel. These arguments will not be repeated 

here. Let us note, however, that for the most part our simulation studies have 

neglected phase fluctuations and have concentrated on the effects of amplitude 

fluctuations alone. The important question in a choice of fading channel model 

is then how closely the Rayleigh-Rice and lognormal distributions resemble each 

2 

other. At least for the relatively small values of indicated in Table 3-l> 

the lognormal distribution is closely approximated by the Rayleigh-Rice distri- 

2 

bution for appropriate choice of a and <5* . This is particularly true in the 
region of small values of r relative to the rms value (r rms = \[z{r 2 } ) which 
essentially determines the perfoi'tnance with coding. In particular, consider the 
exceedence probability given by Q,-,(r) = 1 - F D (r) where F (r) is the cumulative 

A A A 

distribution function (c.d.f.') corresponding to either the lognormal or 
Rayleigh-Rice distribution. For the former we have 

V r) - 1 -5 r ^o (3.17) 

= 1 ; r < 0 


with '<[j (x) the unit Gaussian c.d.f., while for the Rayleigh-Rice distribution we 
have 


^(r) = Q(a, r/<r) 5 

' = i ; 


r £ 0 
r < 0 


(3.18) 
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* L is depth of penetration into Venusian atmosphere 

** B , is 3 db single-side bandwidth of the log- amplitude spectra 


Data based upon homogeneous turbulence below 55 km with structure constant 

, -6 - 1 / 3 

C n = 2.024 x 10 m (inferred from Venera 4) 
and wind shear velocity 

Vj = 50 m/sec and unifora (Venera 4 data) 

Table 3-1 

Summary of Log-Amplitude and Phase Fluctuation Results 
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■where 


CO 

Q(or, p) = | x exp 

P 



x 2 + 


2 


I (a x) dx 
o 


is Marcum’s Q-function (cf. [29]). The quantity Q^(r) is plotted in Figure 3-1 
2 

for selected < 5 ^ and a. It is clear that there is a close resemblence for 

appropriate parameter choices. We have, not attempted to determine an optimum 

choice of parameter values in any sense. Instead we have chosen to evaluate code 

* 2 

performance parametrically as a function of both a and . In particular, 

for fixed X we have equated mean-square values for both the Rayleigh-Rice and 

the lognormal distributions to obtain an empirical relationship between a and 
2 

• Equating mean-square values we obtain 


2 <T 2 + * 2 = exp ^2 


} 


( 3 . 20 ) 


2 2 2 

so that for fixed choice of and X we can solve for CT' ( or equivalently 


(Tl ) with the result 

cl 


o; 2 = 2 CT 2 = exp ^2 (T^} - ^ 


( 3 . 21 ) 


;(i -X 2 ) + 2 <r x 2 


where the approximation is justified due to the relatively small values of < 5 ^ 
as indicated in Table 3 - 1 . It follows that 


a 2 = (§:) = 


2X‘ 


(1-X 2 ) + 2 <r* 


( 3 . 22 ) 


2 

which is the desired empirical relationship. For the small values of < 5 ^ considered 

here the choice X = 1 appears to provide the best match. For example with 

= 0.056 and X = 1 we have from (3-22) a~ 4.2 while ccczz 7.5 for ~ 0.018. 

In either case, the match between the Rayleigh-Rice and lognormal distributions as 

2 

illustrated in Figure 3-1 is quite good and, improves with smaller values of . 



Figure 3-1 

Comparison of Lognormal and Rayleigh-Rice 
Amplitude Distribution 
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Indeed., as argued in [28] it is only for small values of (5 ^ 2 that the lognormal 
distribution for the scattered signal amplitude can be considered valid. 

Although the case X = 1 appears to provide a best fit to the lognormal 
distribution it remains of interest to consider performance over a wide range of 
X values. Due primarily to limitations on computer processing time, however, 
we have been forced to limit our consideration of X values. As a result, we 
have chosen to concentrate on the two cases X = 0 and X = 1* This should 
essentially bound the range of anticipated channel conditions. 

In previous work [2^ , [2^ we choose to restrict attention to channel 
frequency dispersion functions of the form 





2 * B 2 + f 2 
o 


(3.23) 


where B is the coherence bandwidth of the channel in Hz. As a result the auto- 
0 

correlation function of the a(t) process in (3-9) is given by 

p. aa (f ) = E + tr) a*(t)} (3. 24) 

2 ' 

■ 6 

A relationship between B and was derived in [28] with the result 


B o = 'S'* B X 


(3-25) 


p 

With and B-^ specified as in Table 3-1 it is a simple matter to tabulate 

2 

appropriate values of the two important parameters CT* and B . This is provided 

^ Ql O 

in Table 3-2. . 
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55 

0.056 

0.436 

0.146 

1.118 

0.112 

30 

0.018 

0.59 

0.112 

1.037 

0.036 

10 

0.0025 

1.02 

0.071 

1.005 

0.005 

5 

0.007 

1.45 

0.054 

1.001 

0.001 

1 

4 x 10-5 

3.23 

0.029 

1.000 

- 


L is depth of penetration into Venusian atmosphere 
Data taken from Table 3 - 1 • 


Table 3-2 

Summary of Fading Channel Model Parameters 

It is to be recalled that in [25] it was assumed for simulation purposes 
that the process a(t) varies slowly relative to an elementary signaling interval 

of T seconds duration so that it can be considered constant over any such interval 

s • •• . . ■ 

but allowed to vary from interval- to- interval. It was remarked that this is a 
realistic assumption for B T << 1 . From Table 3-2 it is seen that this restricts 
us to signaling rates in excess of approximately one channel symbol/second which 
is felt quite reasonable for planetary entry missions. We will assume that some 
mechanism is available to allow perfect phase tracking of the signal component of 
the received signal. Admittedly this is somewhat unrealistic and later in this 
report we will describe in some detail the effects of imperfect phase tracking. 

For the time being, however, let us assume that a perfect carrier phase reference 
is available. The receiver output, after suitable normalization, is then the 
sequence. 
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r. = x. 
1 1 


2E 

C 

N 


+ a. 

l 


+ n i ’ 


i = 1, 2 ,... 


(3.2 6) 


which replaces that given previously in (2.6). Here a^ represents the value of the 

tYl 

a(t) process throughout the i signaling interval and can be described by the first 
order regression 

a i+1 = ^ a i + W i ’ i = 1, 2,... (3*27) 


The sequence ^ w/^ is an i. i.d. sequence of zero mean complex Gaussian random 
variates such that 

E ^ w.w.^ = 0 (3.28a) 


and 




13 


(3.28b) 


•2«B T 
o s 


and T the duration of an 
s 


with the Kronecker delta function, = e 

elementary signaling interval. This scheme has been described previously in [_2^ 

and need not.be repeated. Simulation studies have been performed parametrically as 

a function of B Q T g and allows conclusions to be drawn as a function of signaling 

rate f = l/T . 
s s 


Finally, it is to be noted, that one of the major concerns in the simulation 
results to be described in the following sections is to establish the efficacy of 
interleaving in combatting the channel memory represented by the time-correlated 
fading.' For short constraint length codes the quantity of interest is the amount 
of interleaving required to achieve a given error probability for specified E^/l^ 
as a function of the parameters describing the channel model, the code and the 
receiver operation. Of particular importance will be the interleaving requirements 
as a function of the channel coherence bandwidth B q with other parameters held 
fixed. For longer constraint length codes we are interested in the effectiveness 
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of interleaving in reducing the computational and/or storage requirements of 
sequential decoders in the presence of fading. While there are many interleaving 
techniques possible, we shall employ a particularly simple approach until such 
time as a superior approach becomes apparent. We shall call this approach block 
interleaving and is described as follows: The encoded serial bit stream is first 

blocked into blocks of length L bits where L = J^n for some Jl = 1, 2, ... and 
n is related to the code rate R = l/n. AnlxL block interleaver then transmits 
each of the L successive bits in a block separated by I - 1 bits from I - 1 other 
blocks. The storage requirements are then N = I x I bits. 

4.0 Simulation Results for Viterbi Decoding 

In this section we will describe some of the simulation results obtained to 
date for short constraint length convolutional codes in conjunction with Viterbi 
maximum likelihood decoding. The channel model will be as described in the pre- 
ceding section. The codes considered in this study are listed in Table 4-1. In 


R bits/channel use 

Constraint Length K 

, 


1/2 

3, 7, 10 

3/3 

3, 6 


Table 4-1 

Short Constraint Length Codes Considered 
in Simulation. Study 


each case we utilize the Odenwalder/Larsen codes (cf. 00, [9]) as discussed pre- 
viously in Section 2.0. The K = 3 code for both rate l/2 and l/3 is included since 

its performance is relatively easy to simulate. In cases where Y = 0 we have 

2 - 

simulated only the case (P =1 since Table 3-2 indicates the actual value will 

cL 

be very close to unity for a wide range of mission parameters. For Y = on the 

2 

other hand, we have considered only the value 0 Tn = 0.1 as this appears repre- 

cl " 

sentative of worst case values appearing in Table 3-2. . 
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Let us first consider the effects of <5^ and/or ^ on the bit error 

s. 

probability when there is little or no memory on the channel. This is obtained 

* 

by setting B q T s equal to some relatively large value, say B q T s = 10.0. The 

• \/ 2 

resulting performance then represents an ultimate limit for fixed 6 and 
which can be approached with sufficiently large interleaving. The simulated per- 
formance, with infinitely fine quantization is illustrated in Figures 4-1 and 4-2 
for X - 0 and = 1.0 and in Figures 4-3 and 4-4 for X = 1-0, = 0.1 

9. o. 

with selected rate l/2 and 1/3 codes. These results should be compared with 
Figures 2-5 through 2-8 for performance over the AWGN channel. It should be clear 
from Figures 4-3 and 4-4 that, for ^ = 1.0 and =0.1, performance close to 

that obtained in the absence of fading can be achieved provided the channel is 
memory less. For smaller values of <5" s , the performance approaches the upper 

3 , 

bound for the AWGN channel even more closely. The case V - 1*0 and <S ^ = 0.1 

cL 

then, according to the discussion of the previous section, provides a useful 
measure of representative performance when the scale of amplitude fading repre- 

2 2 

sented by is small. Although large values of <T^ have not been predicted 

by propagation studies, the case ^ = 0, (S = 1-0 nevertheless provides a 

a 

2 

useful measure of performance for large (£• and in this sense provides a worst 
case performance bound. 

In Figures 4-5 through 4-8 we illustrate the expected performance as a 
function of the dimensionless quantity B q T s when there is no interleaving and 
infinitely fine receiver quantization is employed. It is apparent that severe 
degradation in performance results with increasing channel memory represented by 
small values of B q T s . In particular some degree of interleaving must be provided 
if the signaling rate is such that B T ^0.1. Referring to Table 3-2, for a 

penetration of 55 km into the Venusian atmosphere this would imply a maximum 

— 

Care must be taken in the interpretation of the large B 0 T S results. As stated, they 
provide a useful measure of the ultimate performance which can be approached with 
sufficiently large interleaving. The actual performance for large B 0 T S (i.e., fast 
fading) must take into account the decorrelation loss in the detectors. We have not 
pursued this question. > 


R-'fe, %\*\o.c 
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Figure 4-1 

Simulated Performance of Selected Rate l/2 

Codes with B T - 10.0, Y = 0 and (J' 2 = 1.0 
os a 
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Figure 4-4 


Simulated Performance of Selected Rate l/3 
Codes with B T = 10. 0,# = 1.0 and CT 2 =0.1 
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Figure 4-5 

Simulated Performance of K = 7, R = l/2 Code for 
Different Values of B Q T s without Interleaving 

and with V = 0, = i.o 
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signaling rate of approximately 1 channel symbol per second if interleaving is to 

t 

be avoided. It should be noted that at higher altitudes 0 ^- is smaller and the 
effects of channel memory are not as severe. This is offset, however, by the corre- 
sponding decrease in B . At any rate, the conclusions to be drawn are that some 
degree of interleaving is required if reasonable signaling rates are to be achieved. 

4.1 Interleaving Considerations 

A number of simulations have been performed to establish the efficacy of 
interleaving in reducing the channel memory and thereby avoiding the severe 
degradations in error probability performance for small B q T s demonstrated in the 
preceding section. While considerably more work is required in this area, some 
tentative conclusions can be drawn on the basis of simulation results obtained as 
part of this study. 

An Illustration of typical improvement in bit error probability to be 

realized with the use of a simple block interleaver structure is provided in 

Figures 4-9 and 4-10 for a K = 7, R = l/2 code with B T = 0.001. In all cases 

we have employed an I x L block interleaver as described in Section 3.0 with 

L = 2. n for some n = 1, 2,... and n is related to the code rate by R = l/n. We 

have found it convenient in simulating interleaver performance for convolution 

codes of different rates to parameterize results in terms of the two quantities 

I and $ . As a result we have consistently labeled curves in terms of the literal. 

product 1x2. Note that no interleaving is implied in the case 1 x 2. for any 

i value of 2 . It is clear from Figures 4-9 and 4-10 that, at least for this code, 

an interleaver of size in excess of 30 x 30 is required to approach within 1 db 

of the limiting performance for B q T s in the vicinity of 0.001. Similar conclusions 

can be drawn from Figures 4-11 and 4-12 for a K = 6, R = l/3 code and in Figure 4-13 

for a K = 3> R s l/2 code again with B T = 0.001 in each case. For larger values 

os 
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Effects of Block Interleaving for K = 7, R = l/2 
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Figure k-12 

Effects of Block Interleaving for K = 6 

Code with B T = 0.001 and X = 1.0. 

os 7 a 
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Figure 4-13 

Effects of Block Interleaving for K = 3 > R = l/2 

Code with B T = 0.001 and X = 0, = 1.0 
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of B T , of course, the interleaving requirements are less severe. This is 

OS 

illustrated in Figures 4-l4 and 4-15 for a K = 6, R = l/3 code with B T =0.01 

'os 

and B T =0.1 respectively. The effects of interleaving have shown a surprising 
insensitivity to both code rate and constraint length depending instead, only upon 
the value of B o T g . It is somewhat surprising that such small amounts of inter- 
leaving can be so effective in combatting the memory effects of the channel. 

This behavior is easily explained by reference to Figure 4-l6 which illustrates 
in some detail the approach of the bit error probability performance to the 

limiting case of zero channel memory as a function of B T for a K = 3 > R = l/2 

o s 

code. Observe that for B T = 0.001 initially, a 30 x 30 interleaver results in 

o s 

an "effective" time -bandwidth product of 0.03 which according to Figure 4-l6 
should provide performance within 1 db of the limiting case. In general an inter- 
leaver of size I x £ results in an "effective" time -bandwidth product of 

B T x I since successive symbols are then separated by I channel symbols reducing 
o s 

the memory or equivalently increasing the coherence bandwidth by the factor I. 

This is illustrated in Figure ,4-17 for a K = 7 , R = l/2 code with an initial 
B q T s = 0.001. it is demonstrated that the performance obtained with interleaving 
compares favorably with that obtained with the corresponding "effective" value 
of B Q T g . As a rule of thumb it would appear that if performance within a few 
tenths of a db of the limiting performance is to be obtained, the "effective" 

B q T s should be in the vicinity of 0.1. Table 4-2 provides a summary of the inter- 
leaving requirements as a function of B o T g . 
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Figure 4-15 

Effects of Block Interleaving for K =6, R=l/3 
Code with B T =0.1 and Y -0, 0T^=1.0 
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B T 
o s 

Required I x X 



.1 

.01 

.001 

lxl 
10 x 10 
100 x 100 


Table 4-2 

Summary of Interleaving Requirements as a 

Function of B T to Obtain 
os 

Performance within a Few 
Tenths of a db of 
Limiting Performance 

Finally it is to be noted that in all of the preceding work we have con- 
sidered only square interleaver structures. It is to be expected that performance 
should be an insensitive function of the dimension $ and hence economies can be 
realized by reducing this dimension considerably. Figure 4-l8 illustrates typical 
results for aK = 3, R = l/2 code utilizing non- square interleaving. Obviously 
there is a finite limit to how small $ can be made without seriously degrading 
performance. This is a subject of continuing study. 

I 

4.2 Quantization Effects : 

In all the simulations described so far infinitely fine quantization was 
employed. It remains to determine the effects of receiver quantization. The 
effects of receiver quantization are in clear evidence in Figure 4-19 and Figures 
4-20 and 4-21 for a K = 6, R = l/3 code and B Q T g = 10.0 and B Q T g = 0.01 respectively. 
In each case no interleaving is employed. It appears, as in the case of the AWGN 
channel, that a degradation of only a few tenths of a db result with Q = 8 level 
quantization. We do not expect the previously drawn conclusions on interleaving 
requirements to be altered if 8 level receiver quantization is employed. 
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4.3 Comparison of Results With Lognormal Amplitude Distribution 

While there is some theoretical justification of the lognormal amplitude 

distribution there is also reason to believe that the Rayleigh-Rice amplitude 

distribution provides a more realistic model of actual fading channel behavior. 

It is of some interest to compare simulation results obtained under both assumptions. 

This is illustrated in Figure 4-22 for a K = 3, R=l/2 code and B q T s = 0.001. For 

the lognormal distribution we take = 0.056 and = 0.018 with B-^ related 

to B through (3.25). The parameters of the Rayleigh-Rice distribution have been 
o 

chosen as X = 1.0 and taken from Table 3-2. The simulated bit error proba- 

a. 

bility performance compares favorably with the correspondence improving with smaller 
2 

values of . 


4.4 Imperfect Phase Tracking 

The simulation results reported in preceding sections have all assumed perfect 
phase tracking. As a result, the sequence of matched filter outputs can be described 
by (3-26) and the sequel. In reality, since the phase is not known at the receiver, 
quadrature matched filters must be employed. During the i signaling interval 
the inphase and quadrature matched filter outputs can be respectively represented as 




cos 9. + n T . : 
i I,i 



2 ,. 


(4 . la) 


and 




+ a. 

l 


sin 9. + n. . 
i 



2 ,. 


C'+.lh) 


Here $ and are mutually independent i.i.d 

sequences and 


. zero-mean unit Gaussian 


Bit Error Probability P f 
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Figure h-22 

Comparison of Results for Rayleigh-Rice and Lognormal 
' Distribution with K = 3j H = l/2 Code 
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© i = arg 


[' * <1 


(4.2) 


= tan 


-1 


* sln ^ * a Q„i 

cos 4) + a I ^ i J 5 


i =1, 2,... 


*th 

represents the value of the true phase process Q(t) throughout the i signaling 

interval. The quantities a . and a. . are respectively the inphase and quadrature 

I,i Q,i 

components of the complex quantity a^ described previously, i.e.. 


a i =a l,i + J a Q,i 5 


i = 1, 2, ... 


(4.3) 


as in |25| , It proves convenient to represent the quadrature matched filter outputs 
as 


R i “ r l,i + * r Q,i 


(‘*.4) 


/ 2E , , 38. 

= It i ViT r+ a il e 1 + H i ; 

O 


i = 1, 2, . .". 


where now 

N. = n . + j n . ; i = 1, 2,... (4.5) 

X ± j X X - 

We have then made use of a phase estimator originally proposed hy Jacobs h and 
easily shown to be related to that described by Arnstein [30] . In particular, we 
develop the estimate 


i-1 


i-±i“M Z ®k 2 f 5 


i = 1, 2, . . . 


(*. 6 ) 


. k=i-N 


th 

of the phase ©^ during the i cn signaling interval based upon N past observations of 
the matched filter output sequence represented by (4.4). The 180° phase ambiguity 
represented by the + sign in (4.6) can easily be removed by periodic insertion of 
a known signal sequence and subsequent tracking of the phase during successive 
signaling intervals. An analysis of the estimation accuracy to be achieved with 


this estimator is provided in [ 31 ] where it is shown that for large SNR the 
performance approaches that of a first order phase-locked loop (PLL) with 
tracking bandwidth b = l/NT . The advantage of this approach is that it does 
not require carrier power. If the phase were truly constant for all time it 
would make sense to utilize as large a value of N in (4.6) as possible. The 
assumption of constant phase over the N preceding signaling intervals is 
approximately satisfied if N < 1 /B q T s where the dimensionless quantity B q T s has 
been defined previously. In our simulations we will assume that 


N = 


a 

B T 
o s 


(4 


when 0 < a < 1 is to be specified. The question remains as to how a is to be 
chosen as a function of the fading channel parameters. Figure 4-23 illustrates 
a typical simulation result for a K = 3, R = l/2 code with X = 0, =1.0 and 

B T = 0.015. With a choice Y = 0, the phase varies so rapidly that nothing is 
to be gained by employing N in excess of 3 (<2 = 0.05 in this case). In Figure 
4-24 we illustrate the performance of the phase tracker for a K = 3, R = l/3 code 

O 1 

now with Y = 1.0, (P = 0.1 and B T = 0.001. Here the phase varies so slowly 
that performance within a few tenths of a db of that afforded by perfect phase 
tracking can be achieved with N = 100 (a = 0.1). Additional results are illu- 
strated in Figures 4-25 and 4-26. It would appear that for Y = 1.0, (P^ = 0.1, 
which we feel are typical channel parameters for a Venusian atmosphere, a value a 
in the range 0.25 ^ &. £ 0.1 should provide performance within a few tenths of a 
db of that obtainable with perfect phase tracking. 

5.0 Simulation Results for Sequential Decoding , 

In this section we will describe some of the simulation results obtained 
for both Fano and ZJ decoding of long constraint length convolutional codes. We 
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Figure 4-23 

Effects of Imperfect Phase Tracking on K = 3, R = l/2 
Code with X = 0, CP- 2 = 1.0 and B T g = 0.015 
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will attempt to parallel the discussion of the Viterbi decoding results presented 
in the preceding section to the extent possible. It is to be emphasized at the 
outset that considerably less simulation results have been obtained for sequential 
decoding vis-a-vis the results for Viterbi decoding. This is due chiefly to the 
time consuming and expensive nature of the computer processing required in order 
to obtain statistically significant data. As a consequence, the results and con- 
clusions reported here should be considered tenative based as they are on a 
limited number of simulation results. It is expected, however, that all sub- 
stantive issues will be resolved in the course of work presently in progress. 

In all simulations of sequential decoder performance we have made use of 
the Massey quick-look code described previously in Section 2 with constraint length 
K = 32 and rate R = l/2. Each frame or block consisted of 250 information bits 
to which 31 tail zeros -were added. The receiver employed uniform quantization 

4 

with Q = 8 levels. All simulation results are reported on the basis of 10 

successive transmissions of the 250 bit message. A computational cutoff of 

4 

5 x 10 computations was imposed for both the Fano and ZJ decoder. The ZJ 

• 3 

algorithm employed a decode table consisting of a maximum of 5 x 10 entries 
although only a small percentage of this number were ever required. Both the 
Fano and ZJ decoders employed a bias \ = l/2 (i.e., equal to the rate) while 

the Fano decoder employed a threshold spacing T q = 4 bits. We have found it 
convenient to present simulation results for fixed values of r/R q . Here R q 
refers to the computational cutoff rate or critical rate in the absence of fading, 
i.e., in the AWGN channel as determined by (2.13) and the sequel. It is recog- 
nized that the actual critical rate is effectively lowered in the presence of 
fading although no attempts have been made to explicitly compute this quantity 
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as a function of the underlying channel parameters. In what follows then we will 

continue to use the normalized rate R/R q referred to the critical rate of the 

original AWGN channel as a convenient parameter in describing the behavior of 

sequential decoders in the planetary entry channel. 

In Figures 5-1 and 5-2 we illustrate the empirically determined computational 

distribution of the Fano decoder as a function of the dimensionless quantity B T 

os 

for R/r = 0.8 and 0.9 respectively. Here the channel parameters were chosen as 
o 

V= 1.0 and CP 2 =0.1. For the case R/R = 0.99, the probability of exceeding the 
a ' o 

4 

imposed computational cutoff of 5 x 10 computations per block was close to unity 
even for B T = 10.0 so that results are not plotted. Similarly, for the choice 

OS . • .. • # _ 

2 

of channel parameters V = 0.0 and (J^ =1.0 the quit probability was con- 
sistently found to be close to unity for r/R q ^ 0*8* As a result, this worst 
case has been eliminated from further consideration. On Figures 5-1 and 5-2 we 

have indicated the performance in AWGN for comparison purposes. The case B T =10.0 

o s 

represents, for all practical purposes, the limiting case of zero channel memory 
which can be approached with sufficiently large interleaving. It is clear that 
even moderate amounts of channel fading cam result in severe degradation of the 
computational performance of the Fano decoder. Furthermore, even if large amounts 
of interleaving is provided the best that one could expect is to approach the 
limiting case of zero channel memory represented by the case B q T s = 10.0. For 
r/R q close to unity, this performance can in turn be orders of magnitude worse 
than the computational requirements in the AWGN channel. Similar comments can be 
made for the ZJ decoder on the basis of the simulation results obtained 
to date. In the absence of fading, prudent system design would dictate 

choosing r/R q close to unity. From the limited simulation results reported here 
it must then be concluded that, in the presence of even moderate amounts of channel 
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Figure 5-1 

Empirical Probability Distribution of Computational Count 
for Fano Decoder with V = 1.0, (T? ~ 0.1 and r/r = 0.8 
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fading, interleaving alone will not prove sufficient to insure computational 
requirements comparable to that of the AWGN channel. What is required is clearly 
a combination of interleaving and a simultaneous reduction in the operating r/R q 
ratio. While it is felt that the required reduction in R /r can be explicitly 
calculated, no attempt was made to do so as part of this study. It is felt that 
this is a topic worthy of additional investigation. 

5.1 Interleaving Considerations 

A limited number of simulations have been performed to investigate the 

efficacy of simple block interleaving in combatting the channel memory and thereby 

reducing somewhat the computational requirements of sequential decoders for small 

values of B T . Typical results for a Fano decoder are illustrated in Figures 5-3 
os / • 

and 5-^ for R/R q = 0*8 and 0.9 respectively with B q T s = 0.01 in either case. 

Here it would appear that an interleaver of size 50 x 50 will suffice to provide 
performance comparable to the memory less channel if R/R q = 0.8 while for 
R'/R = 0-9 an interleaver of size 100 x 100 is required. Unfortunately time and 
computer processing costs have precluded obtaining additional simulation results. 

It is recommended that future investigations concerning the efficacy of inter- 
leaving consider the more meaningful range of r/R q values such that r/R q 0.8. 

6.0 Summary and Conclusions 

An attempt has been made to describe an approach to the determination and 
parameterization of the performance of convolutional codes on fading channels 
typical of planetary entry missions. For short constraint length codes and Viterbi 
decoding we have considered the effects of amplitude fading alone under the assumption 
a perfect carrier phase reference is available and have shown the severe degradation 
which results unless some form of interleaving is utilized to combat the channel 
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memory. The efficacy of simple block interleaving was investigated and shown to 
be quite effective in reducing the bit error probability P^. The effects of 
imperfect carrier phase tracking was then considered and a simple phase estimator 
proposed which has shown to result in performance comparable to that of perfect 
phase tracking for reasonable E^/n and provided the parameters of the phase 
tracker were chosen appropriately. Other than some additional work on the effects 
of phase tracking errors, it is felt that the treatment of short constraint length 
codes and Viterbi decoding is complete. 

For longer constraint length codes and sequential decoding much more work 
remains. In particular, additional work should be done to establish the efficacy- 
of block interleaving than has been reported here. The penality in operating 
r/R q to maintain computational and/or storage requirements comparable to that for 
the AWGN channel should be investigated in detail. We have not studied the 
effects of imperfect phase tracking to any great extent and this subject is 
certainly worthy of detailed investigation. Finally, it appears that with 
suitable implementation, the ZJ decoder is capable of .providing near real-time 
decoding of long constraint convolutional codes. This could be of considerable 
utility in relay links where decoding is provided onboard the spacecraft. Much 
work remains, however, before a complete understanding is available concerning 
the nature of the tradeoff between the block error probability P^, the quit error 
probability P n and the computational and storage requirements. These topics are 
all subjects of continuing investigation. 
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